Submitted to Ap.J. Letters, June 16, 2000 
Preprint typeset using LAT^X style emulateapj v. 26/01/00 



O 
O 



in 

(N 
> 

in 

(N 

\o ■ 
o : 
o . 
o ■ 

6 ■ 
a: 

C/3 . 



13 



THE AMPLITUDE OF THE KILOHERTZ QUASI-PERIODIC OSCILLATIONS IN 4U 1728-34, 4U 1608-52, 
AND AQL X-l, AS A FUNCTION OF X-RAY INTENSITY 

Mariano Mendez 1,2 , Michiel van der Klis, Eric C. Ford 

Astronomical Institute "Anton Pannekoek", University of Amsterdam and Center for High-Energy Astrophysics, Kruislaan 

403, NL-1098 SJ Amsterdam, the Netherlands 

Submitted to Ap.J. Letters, June 16, 2000 
ABSTRACT 

We study the kilohertz quasi-periodic oscillations (kHz QPOs) in the low-mass X-ray binaries 4U 1728-34, 4U 
1608-52, and Aql X-l. Each source traces out a set of nearly parallel tracks in a frequency vs. X-ray count rate 
diagram. We find that between two of these tracks, for similar QPO frequency, the 2 — 60 keV source count rate 
can differ by up to a factor of ~ 4, whereas at the same time the rms amplitude of the kHz QPOs is only a factor 
of ~ 1.1 different. We also find that, for 4U 1608-52 and Aql X-l, the rms spectrum of the kHz QPOs does 
not depend upon which track the source occupies in the frequency vs. X-ray count rate diagram. Our results for 
4U 1728-34, 4U 1608-52, and Aql X-l are inconsistent with simple "extra source of X-rays" scenarios for the 
parallel tracks, such as those in which the properties of the kHz QPOs are only determined by the mass accretion 
rate through the disk, whereas X-ray count rate also depends upon other sources of energy that do not affect the 
QPOs. 



Subject headings: accretion, accretion disks - 
Aql X-l) — X-rays: stars 

1. INTRODUCTION 

It is thought that the kilohertz quasi-periodic oscillations 
(kHz QPOs) in low-mass X-ray binaries reflect the motion of 
matter in orbit at some preferred radius in the accretion disk 
around the neutron star (Miller, Lamb, & Psaltis 1998; Stella & 
Vietri 1999; Osherovich & Titarchuk 1999; Cui 2000; Cam- 
pana 2000). In most sources, two simultaneous kHz QPOs 
have been detected, with frequencies between 300 and 1300 
Hz. In a given source, the frequency of the QPOs can shift 
by a few hundred Hz, apparently as a function of M, the rate 
at which mass is accreted onto the neutron star (Miller et al. 
1998). Yet, it is well known that the relation between the kHz 
QPO frequencies and X-ray count rate (or X-ray flux), which is 
assumed to be a good indicator of Al, is complex (Ford et al. 
1997a; Zhang et al. 1998; Mendez et al. 1999). In a frequency 
vs. intensity {v — Ix) diagram a source displays a set of almost 
parallel tracks. On time scales of a day or less the source moves 
along one of these tracks with QPO frequency and X-ray count 
rate positively correlated, whereas on long time scales, in ob- 
servations separated by a few days the source occupies different 
tracks in this diagram (see, e.g., Fig. 2 in Mendez et al. 1999). 

On time scales of a few days, QPO frequency correlates 
much better with the X-ray spectral properties of the source 
than with count rate. In 4U 0614+09 it was found that the fre- 
quencies of the kHz QPOs correlate better to the flux of the soft 
("blackbody") component in the energy spectrum than to the to- 
tal X-ray flux (Ford et al. 1997b), whereas a good correlation 
of QPO frequency with photon index of the power-law com- 
ponent in the energy spectrum has been found in 4U 0614+09 
and 4U 1608-52 (Kaaret et al. 1998). The frequencies of the 
kHz QPOs also correlate much better with the position of the 
source in an X-ray color-color diagram than with count rate 
(Wijnands et al. 1997, 1998b; Jonker et al. 1998; Mendez et 
al. 1999; Mendez & van der Klis 1999; Kaaret et al. 1999a; 
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Mendez 2000; Reig et al. 2000; van Straaten et al. 2000; Di 
Salvo et al. 2001). This correlation extends to other timing 
properties as well. The power spectra of these sources often 
show a broad-band noise component that is roughly flat below 
a break frequency at v\> ~ 1 — 10 Hz and decreases above v-^,, 
and a low-frequency QPO at ~ 10 — 50 Hz. Studies with EX- 
OSAT and Ginga already demonstrated that the rms amplitude 
and break frequency of the broad band noise component, and 
the frequency and rms amplitude of the low-frequency QPO 
depend monotonically upon the position of the source on the 
color-color diagram, and correlate much better to colors than to 
count rate (see van der Klis 1995, for a review). Studies with 
the Rossi X-ray Timing Explorer (RXTE) confirmed this across 
a much wider range of sources (Wijnands et al. 1997, 1998a,b 
Jonker et al. 1998; Markwardt, Strohmayer, & Swank 1999 
Boirin et al. 1999; Jonker et al. 2000; Di Salvo et al. 2001 
Bloser et al. 2000), and also showed that both v\, and the fre- 
quency of the low-frequency QPO are strongly correlated to the 
frequency of the kHz QPOs (Stella & Vietri 1998; Ford & van 
der Klis 1998; Psaltis, Belloni, & van der Klis 1999; Di Salvo 
etal. 200 1 ) and to each other ( Wijnands & van der Klis 1999). 

These correlations have reinforced the idea (Hasinger & van 
der Klis 1989) that there is one single parameter, inferred to 
be Al, which governs all the timing and spectral properties of 
these sources; X-ray flux, which can change more or less in- 
dependently of the other parameters, is the exception (see also 
van der Klis et al. 1990; Mendez et al. 1999). 

In principle, the disparate behavior of X-ray flux and M 
could be understood if mass flows onto the neutron star through 
two independent channels (e.g., Ghosh & Lamb 1979; Former, 
Lamb, & Miller 1989). Part of the mass accretion onto the star 
comes from the disk flow, Md, whereas the rest of the mass 
flux onto the star comes from a radial flow, Mr. The timing 
and spectral properties are governed by Md, whereas the X-ray 
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FIG. 1 . — The amplitude of the kHz QPOs as a function of their frequency 
for (a) 4U 1728-34, (b) 4U 1608-52, and (c) Aql X-l. For 4U 1728-34 and 
4U 1608-52 the filled and open circles show the relation for the lower and 
the upper kHz QPO, respectively, and open squares indicate measurements in 
which we only detect a single QPO. For Aql X-l the open squares show the 
relation for the only kHz QPO so far detected. The rms amplitudes are for the 
full PCA energy band. 

luminosity is proportional to M D + Mr (Wijnands et al. 1996; 
Kaaret et al. 1998; Mendez et al. 1999). In this scenario, 
the frequencies of the kHz QPOs are linked to Md, but because 
Md and Mr can vary more or less independently, and both con- 
tribute to the total luminosity (Ld + £r)> the source traces out 
separate tracks in the v — Ix diagram. (The argument is more 
general, and is still valid if Lr is produced by another source 
of energy, not necessarily related to a mass inflow.) None of the 
kHz QPO models so far discussed explicitly predicted the exis- 
tence of the parallel track phenomenon, and there are therefore 
no specific predictions for how QPO amplitude varies across 
tracks. For this reason in this paper we use as a reference the 
simplest possible scenario, namely, one in which the additional 
source of energy that produces the extra luminosity producing 
the parallel tracks does not affect the QPO amplitudes at all. 
In such a scenario, at a constant kHz QPO frequency (i.e., con- 
stant M D ), the rms amplitude of the QPOs should decrease with 
count rate, as the source moves from one track to another in the 
v — Ix diagram. 

In this paper we use RXTE data to study the rms fractional 
amplitude of the kHz QPOs on different tracks in the v — Ix 
diagram, that is, as a function of count rate for a fixed QPO 



frequency. In §2 we describe the observations and results, and 
in §3 we discuss the significance of our findings within the sce- 
nario previously described. 

2. OBSERVATIONS AND RESULTS 

We used observations from the Proportional Counter Array 
(PCA) on board RXTE of the kHz QPO sources 4U 1608- 
52, 4U 1728-34, and Aql X-l. For 4U 1608-52 the data are 
the same as those in Berger et al. (1996) and Mendez et al. 
(1998a), except that here we do not use the observation of 1996 
March 6 because only photons from a limited energy range 
were recorded. For 4U 1728-34 we analyzed the same data 
as described in Strohmayer et al. (1996, 1997), Ford & van der 
Klis (1998), Mendez & van der Klis (1999), and Di Salvo et 
al. (2001). For Aql X-l we analyzed the same observations as 
used by Zhang et al. (1998); Cui et al. (1998), and Reig et 
al. (2000). We discarded those observations for which fewer 
than five detectors of the PCA were on. This means a loss of 
~ 2 %, ~ 2 %, and ~ 11 % of the observations of 4U 1728-34, 
4U 1608-52, and Aql X-l, respectively. A few X-ray bursts oc- 
curred during these observations, but we excluded those parts 
of the data from the rest of our analysis. 

In all these observations data were recorded using the Stan- 
dard 2 mode, with 16 s time resolution in 129 energy channels, 
and in all cases another mode with a time resolution of at least 
1/4096 s was running in parallel. We used the Standard 2 data to 
produce background subtracted light curves (using Pcabackest 
v2.1b) for the full energy range of the PCA (nominally 2 — 60 
keV). 

For each source we used the high-time resolution data to pro- 
duce power spectra every 64 s, up to a Nyquist frequency of 
2048 Hz, with no energy selection. We searched these power 
spectra for kHz QPOs, at frequencies >250 Hz. To do this 
we first produced a dynamic power spectrum of each observa- 
tion, we identified those power spectra that showed a strong 
QPO, and we visually measured its frequency, Vq. If two QPOs 
were present in a power spectrum, we always picked the one at 
lower frequencies, because it was generally the most conspicu- 
ous one. We fitted those power spectra, in the range v — 100 Hz 
to i>o + 100 Hz, using a function consisting of a constant plus 
one Lorentzian. In a few cases in which the QPOs were not sig- 
nificant enough (less than 3<r, single trial) in the 64 s segments 
we averaged several contiguous 64 s segments, but always less 
than 20, to try and measure the QPO. We did not combine more 
than 20 consecutive power spectra because the QPO frequen- 
cies can vary by a few tens of Hz within a few hundred seconds, 
and these frequency variations result in biased measurements of 
the QPO frequencies. We discarded those data for which this 
procedure did not reveal any significant kHz QPO. We note that 
in this manner we might have discarded data with weak QPOs 
that could have been detected using more data. 

On several occasions we also detected the upper kHz QPO, 
either directly in the dynamical power spectrum, or after we av- 
eraged (or shifted and averaged; see below) the power spectra 
of an observation. In those cases, to increase the significance of 
this QPO and to facilitate measuring it, we aligned the individ- 
ual power spectra according to the frequency of the lower QPO, 
collected them into groups such that the frequency of this QPO 
did not vary by more than 10 — 50 Hz, and calculated the aver- 
age power spectrum for each group (see Mendez et al. 1998b). 
We fitted these power spectra, in a range of 200 Hz around the 
approximate frequency of the upper kHz QPO, using a constant 
plus one Lorentzian. 
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FIG. 2. — The frequency of the lower kHz QPO vs. the 2 — 60 keV count 
rate for (a) 4U 1728-34, (b) 4U 1608-52, and (c) Aql X-l. The count rates 
have been corrected for background, and are for 5 detectors. For 4U 1728-34 
and 4U 1608-52 we only include data where both kHz QPOs were detected 
simultaneously, so we can unambiguously identify the lower kHz peak. For 
Aql X- 1 we plot the only kHz QPO detected in the power spectrum. For 4U 
1608-52 and Aql X-l we indicate 3 tracks on which we calculate the photon 
spectrum of the QPO (see text and Fig. 7). 

In those cases in which we only detected one QPO, we ana- 
lyzed the average power spectrum of that observation to try and 
identify it as either the lower or the upper kHz QPO. If the av- 
erage power spectrum still showed a single QPO, we calculated 
the shifted-and-averaged power spectrum, using the frequency 
of the QPO we had detected as a reference to shift the indi- 
vidual power spectra. If no other QPO was detected after this 
procedure (this was the case for all Aql X-l data; see below) 
we treated that observation separately, so that in the rest of our 
analysis the identification of the kHz QPOs as upper or lower is 
unambiguous. 

In Figure 1 we plot the rms fractional amplitude of the QPOs 
that we detect in each source, as a function of their frequency 
(see also Di Salvo et al. 2001, Fig. 3b). For 4U 1728-34 
and 4U 1608-52, when we detect two kHz QPOs simultane- 
ously, we use filled and open circles for the rms amplitude of 
the lower and the upper kHz QPO, respectively; when we only 



detect a single QPO we use open squares to plot its rms ampli- 
tude. For Aql X-l we use open squares to plot the rms ampli- 
tude of the only kHz QPO so far detected. It is apparent that in 
4U 1728-34 and 4U 1608-52 the amplitude of the upper kHz 
QPO is usually larger than that of the lower kHz QPO, and that 
the amplitude of the upper kHz QPO decreases monotonically 
with frequency, whereas the rms amplitude of lower kHz QPO 
first increases, and then decreases with frequency. When only 
one QPO is detected in 4U 1728-34 and 4U 1608-52, the rms 
amplitude of that QPO (open squares) is consistent with the re- 
lation for the upper kHz QPO in those sources. 

Aql X-l is different from 4U 1608-52 and 4U 1728-34 in 
that it only shows a single kHz QPO in its power spectrum 
(Zhang etal. 1998; Cui et al. 1998; Reig et al. 2000). We ap- 
plied the procedure described above to each of the observations 
of Aql X-l in which this kHz QPO was detected, but in none 
of them could we detect a second QPO. Because in Aql X-l we 
only see a single kHz QPO, it is difficult to tell whether in all 
cases it is the same peak, and if so, whether it is the lower or 
the upper kHz QPO. However, (i) plots of the QPO frequency 
vs. X-ray color (Mendez 2000), or vs. position of the source 
in a color-color diagram (Reig et al. 2000) are consistent with 
one single relation; (ii) the average FWHM of the QPO in Aql 
X-l is ~ 10 Hz (it goes from ~ 2 to ~ 20 Hz), similar to the 
FWHM of the lower kHz QPO (but significantly less than the 
average FWHM of the upper kHz QPO) in 4U 1728-34, 4U 
1608-52, and 4U 1636-53 (Strohmayer et al. 1996; Mendez et 
al. 1998a; Mendez, van der Klis, & van Paradijs 1998c; Kaaret 
et al. 1999b); (iii) the relation of the rms amplitude of the 
QPO vs. the QPO frequency is similar to that of the lower kHz 
QPO in 4U 1608-52 and 4U 1728-34 (see Fig. 1); (iv) the pho- 
ton energy spectrum of the kHz QPO in Aql X-l (see below) 
is consistent with that of the lower kHz QPO in 4U 1608-52, 
and significantly different from that of the upper kHz QPO in 
the same source (Berger et al. 1996; Mendez et al. 1998b). 
We conclude that whenever a kHz QPO was detected in these 
observations of Aql X-l, it was always the same one, and we 
identify it as the lower kHz QPO. In what follows we use all the 
data with the kHz QPO from this source. 

Figure 2 shows, for 4U 1728-34, 4U 1608-52, and Aql X-l, 
the relation between the frequency of the lower kHz QPO, v\ , 
and the 2 — 60 keV background subtracted count rate. As we 
explained above, in a few occasions we had to average several 
consecutive 64 s segments to get a significant measurement of 
the QPOs, and for that reason each point in this figure repre- 
sents between 64 (~ 63% of the points) and 1280 s of data 
(only ~ 8 % represent more than 256 s of data). The figure 
shows, for each source, the well-known parallel tracks in the 
v - I x diagram (e.g., Ford et al. 1997b; Mendez et al. 1999). 

For each source, we measured the rms fractional amplitude 
of the lower and the upper kHz QPO, at constant QPO frequen- 
cies, as a function of count rate. For this purpose we grouped 
the power spectra according to the frequency of the lower QPO 
(Fig. 2), such that v\ did not vary by more than ~ 40 Hz in 
each frequency range. For every source, we averaged the power 
spectra in each frequency range on each of the more or less par- 
allel tracks in Figure 2. For every frequency and count rate in- 
terval, the average power spectrum was calculated by first align- 
ing the individual power spectra according to the frequency of 
the lower QPO displayed in Figure 2. In this way, if there was 
a second QPO in the power spectra, it would become more sig- 
nificant after the alignment, and would be easier to measure its 
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FIG. 3. — a. Rms amplitude of the lower kHz QPO in 4U 1728-34 as a function of the 2 — 60 source count rate, for different lower QPO frequency intervals as 
indicated. The rms amplitudes are for the full PCA energy band. b. Same as Fig. 3a, for the upper kHz QPO in 4U 1728-34. 



properties (cf. Mendez 2001). 

We fitted each of these average power spectra (one per fre- 
quency interval and per track for each source) using a constant 
plus one or two Lorentzians. The errors in the fit parameters 
were calculated using A\ 2 = 1 (la for a single parameter). 
For 4U 1608-52 and 4U 1728-34, when the upper QPO was 
not significant, we fixed its frequency to v\ + Av, with Av de- 
termined from Mendez et al. (1998a, Fig. 3) and Mendez & 
van der Klis (1999, Fig. 1), respectively, and we calculated 
upper limits to its amplitude using A\ 2 = 2.71, corresponding 
to a 95 % confidence level. 

Figures 3 and 4 show the rms amplitude of the lower and the 
upper kHz QPOs as a function of count rate, for different fre- 
quency intervals in 4U 1728-34 and 4U 1608-52, respectively; 
similarly, Figure 5 shows the rms amplitude of the only kHz 
QPO in Aql X-l. Each panel in those figures presents data for 
the indicated frequency interval; each point belongs to an indi- 
vidual track in Figure 2. From these figures it is apparent that, 
for a fixed QPO frequency, the source count rate can increase 
by up to a factor of <~ 4, but the rms fractional amplitude only 
decreases by a factor of ~ 1.1. There are 13 cases for which 
the rms is consistent with being constant, and 9 cases in which 
it decreases with count rate. 

Assuming that a fraction £ of the X-ray count rate Id pro- 
duced by the disk mass inflow Md is modulated, the frac- 
tional rms amplitude of a QPO can be expressed as rms = 
£/d/ (Id + Ir), where /r, is the X-ray count rate produced by 
Mr. If all QPO properties are a function of In only, at constant 
QPO frequency the ratio of the QPO rms amplitudes on tracks 
j and k in the v — J x diagram, rmSj/rms^, has to be equal 
to the inverse ratio of count rates in those tracks, Ix.k/Ix.j- 



Except in the case of the upper kHz QPO in 4U 1608-52 for 
which only two measurements were available, we fitted a lin- 
ear relation to the rms amplitude ratio vs. the inverse count 
rate ratio, rms-Jrmsy = alx,k/lx,j + b, for each QPO in- 
dependently. For the lower kHz QPQ in 4U 1728-34 we get 
a = 0.16 ± 0.09 and b = 0.83 ± 0.07, respectively, xl = 2 for 
27 degrees of freedom. For the upper kHz QPO in 4U 1 728-34, 
a = -0.41 ±0.41 and b = 1. 23 ±0.36, respectively, X 2 u = 0.9 
for 7 degrees of freedom. For the lower kHz QPO in 4U 1608- 
52, a = 0.25±0.06and6 = 0.75±0.04, respectively, X l = 5.1 
for 27 degrees of freedom. For the only kHz QPO in Aql X- 1 , 
a = 0.53±0.09and6 = 0.50±0.07, respectively, xl = 1-1 for 
20 degrees of freedom. (Some of the fits are formally unaccept- 
able which could be due, for instance, to frequency-dependent 
changes of the QPO rms.) Except for the upper kHz QPO in 
4U 1728-34, for which the error bars are too big to draw any 
conclusion, for all the other QPOs both a and b are significantly 
different from the expected values, a = 1 and 6 = 0. 

In Figure 6 we summarize this information. We show the 
ratio of the kHz QPOs rms amplitude vs. the inverse of the 
corresponding count rate ratio, for points from the individual 
panels in Figures 3, 4, and 5 (i.e., for a small frequency range 
of the lower kHz QPO) taken in pairs (we only used data from 
panels that contain 2 or more points in them). The dashed 
line in that figure shows the relation expected in simple "ex- 
tra source of X-rays" scenarios, rmsj/rms^ = /x,k/^x,j- The 
solid line indicates the best fit to the data using a straight line, 
rmsj/rmsk = alx,k/lx,j + b. The best fit parameters are 
a = 0.28 ± 0.04 and'& = 6.73 ± 0.03, respectively (the uncer- 
tainties are the la errors) with a reduced \ 2 of 3 for 92 degrees 
of freedom (which indicates that a single line is not a good fit 



M. Mendez et al. 



5 



10 



9 - 



° 9 - 



.o - 



7 - 



640 - 680 Hz || 680 - 720 Hz 10 - 



K 9.5 P 4 720 - 760 Hz 4/ . 760 - 800 Hz- ^"10 



I* 



800 - 840 Hz 



m 6 



640 - 680 Hz 

:ii ; 

'- + : 


+ -J" 680 - 720 Hz 

^ I 

: ii ; 


1 1 1 1 1 1 1 1 1 1 1 1 

r<!U — fbU riZ 

: I 


1 1 1 1 1 1 1 1 1 1 1 1 

rbU — oUU rlZ 

: i T ; 

T 


i i i i i i i i i i i i i i 

: I : 

: 1 i : 


' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

> R4-0 H7 
T - 



1000 2000 3000 1000 2000 3000 1000 2000 3000 1000 2000 3000 

2 — 60 keV count rate (c/s) 2-60 keV count rate (c/s) 

FIG. 4.— a. Same as Fig. 3a, for the lower kHz QPO in 4U 1608-52. b. Same as Fig. 3a, for the upper kHz QPO in 4U 1608-52. 



: < 700 Hz : 

T 

" i i i i i i i i i i i i i i i i i i i r 


• 740 - 780 Hz : 

' 1 : 

: T ; 

r i i i i i i i i i i i i i i i i i i i i " 


: 780 - 820 Hz : 


: 820 - 860 Hz : 

L - * - 


~ i i i i i i i i i i i i i i i i i i i r 
: 860 - 900 Hz 1 

; ,m ; 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

7 > 900 Hz ~ 
f f 



800 1600 2400 800 1600 2400 
2—60 keV count rate (c/s) 



FIG. 5.— Same as Fig. 3a, for the only kHz QPO in Aql X-l. 

to the data, as might be expected if, for example, there are some 
source to source changes). We also fitted the data to a constant, 
which yields a reduced \ 2 of 4.8 for 93 degrees of freedom. 
This means an F-test probability of 1 x 1CP 10 against the hy- 
pothesis that the fits using a linear function and a constant are 
statistically equivalent. In the case of the linear fit, both a and 
b are significantly different from their expected values. The re- 
duced x 2 is significantly larger (54 for 94 degrees of freedom) 
if we fixed a = 1 and 6 = 0. 

Finally, we measured the rms energy spectrum of the lower 
and the upper kHz QPOs in 4U 1608-52, and the only kHz 
QPO in Aql X- 1 , when these two sources occupied three differ- 
ent tracks in the v — Ix diagram (Fig. 2), at constant lower QPO 
frequency. (Unfortunately, during the 1996 observations of 4U 



1728-34 -the points that trace out the track at the highest count 
rate in Fig. 2a- the high-time resolution data contain no energy 
information, so that it was not possible to do the same analysis 
for this source.) We divided the data into four energy bands, 
2.0 - 4.6 - 8.0 - 13.0 - 21.5 keV, and calculated the power 
spectra for both sources in each band. In Figure 7 we show the 
rms energy spectrum of these QPOs. For 4U 1608-52 we se- 
lected data with QPO frequency between 700 Hz and 800 Hz, 
from tracks 1, 2, and 3 in Figure 2b, at ~ 1300 — 1500 counts 
s~\ - 1600 - 1800 counts s" 1 , and - 1800 - 2400 counts 
s _1 , respectively; for Aql X-l we selected data with QPO fre- 
quency between 820 Hz and 860 Hz, from tracks 1, 2, and 3 in 
Figure 2c, at - 1000 counts s~\ - 1200 - 1400 counts s~\ 
and ~ 1800 — 2200 counts s _1 , respectively. For each source, 
we use different symbols to indicate measurements at each of 
these tracks (see figure caption). Because in 4U 1608-52 the 
upper kHz QPO is less significant than the lower kHz QPO, for 
the upper QPO we measured the rms spectrum using the com- 
bined data of tracks 1, 2, and 3. Figure 7 shows that (i) in both 
sources, the rms energy spectrum of the QPO is independent of 
the track that the source occupies in the v — ix diagram (Fig. 
2; notice that the overall normalization of the rms spectrum de- 
creases as count rate increases, corresponding to the decrease 
of rms amplitude with count rate shown in Fig. 4 and 5); (ii) 
the rms spectrum of the only kHz QPO in Aql X-l is similar 
to that of the lower kHz QPO in 4U 1608-52, but significantly 
different from that of the upper kHz QPO in 4U 1608-52. 

3. DISCUSSION 

For 4U 1728-34, 4U 1608-52, and Aql X-l the frequencies 
of the kHz QPOs trace several, almost parallel, tracks when 
plotted vs. X-ray count rate (Fig. 2). We measured the rms 
fractional amplitude of the kHz QPOs in these three sources at 
different values of the QPO frequency, as a function of count 
rate. We find that, for a fixed QPO frequency, moving from 
one track to another in Figure 2, while the count rate increases 
by up to a factor of <~ 4, the rms fractional amplitude of the 
QPOs only decreases by a factor of ~ 1.1 (Fig. 3, 4, 5, and 
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FIG. 6. — Ratio of the kHz QPO rms amplitude vs. the inverse of the cor- 
responding count rate ratio, for each pair of points from the panels with two 
or more points in Figures 3, 4 , and 5. The lower and upper kHz QPOs in 
4U 1728-34 are plott ed using black and gray filled circles, respectively. The 
lower and upper kHz QPOs in 4U 1608-52 are plotted using black and gray 
open circles, respectively. The kHz QPO in Aql X-l is plotted using black 
open stars. The dashed line shows the expected relation for the "extra source 
of X-rays" scenario with slope 1 and intercep t (see text). The solid line 
is the best fit straight line to the data, with slope 0.28 ± 0.04 and intercept 
0.73 ±0.03. 

6). This result is inconsistent with simple "extra source of X- 
rays" scenarios in which the separate tracks in Figure 2 arise 
from a luminosity component that does not affect the QPO. Fur- 
thermore, we also find that the rms spectrum of the lower kHz 
QPO in 4U 1608-52 and that of the only kHz QPO in Aql X-l 
do not depend upon which track the source occupies in Figure 
2. 

It is well known that the lack of correlation between kHz 
QPO frequencies and X-ray count rate on time scales longer 
than a day (Fig. 2) can coexist with a very good correlation 
between frequency and position in the X-ray color-color dia- 
gram (see §1). It was proposed that this dichotomy could be 
explained in terms of two independent mass accretion rates, a 
disk and a radial inflow; while both of them determine the lu- 
minosity of the source, only the mass flow through the disk 
determines the frequency of the kHz QPOs (Kaaret et al. 1998; 
Mendez et al. 1999). Our results for 4U 1728-34, 4U 1608- 
52, and Aql X-l are inconsistent with version of this scenario 
where the variations in the radial flow do not appreciably affect 
the QPO amplitude: the rms fractional amplitudes of the kHz 
QPOs measured at constant frequency do not decrease rapidly 
enough as count rate increases (Fig. 6). 

The changes in the X-ray intensity may be due to anisotropies 
in the radiation pattern. It is, however, unlikely that the 
radiation that is modulated in the kHz QPO mechanism is 
anisotropic, as kHz QPOs have been observed at roughly simi- 
lar amplitude in sources spanning a large range of inclinations 
(e.g., Barret et al. 1997; Homan & van der Klis 2000). Per- 
haps, there is a large-scale redistribution of some of the emitted 
energy over unobserved spectral bands. Alternatively, the miss- 
ing energy might be used to accelerate a jet. While radio emis- 
sion has been detected from Aql X-l during X-ray outbursts, 
and perhaps from 4U 1728-34, it is unlikely that this explana- 
tion applies to 4U 1608-52, as no radio emission has ever been 
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FIG. 7. — Rms energy spectrum of (a) the lower and upper kHz QPOs in 
4U 1608-52, and (b) the only kHz QPO in Aql X-l. For the lower QPO in 
4U 1608-52 we use filled black, open black, and filled gray circles to indicate 
measurements from tracks 1, 2, and 3 in Figure 2b, respectively. We use gray 
open squares to indicate measurements of the upper kHz QPO in 4U 1608-52. 
In this case we use the combined data from tracks 1, 2, and 3. For the QPO 
in Aql X-l we use filled black, open black, and filled gray triangles to indicate 
measurements from tracks 1,2, and 3 in Figure 2c, respectively. 

detected from this source (Fender & Hendry 2000). 

A similar parallel tracks phenomenon is observed in the en- 
semble of sources with kHz QPOs (van der Klis 1997; Ford et 
al. 2000a, see also Zhang et al. 1997): While for individual 
sources, on a time scale of hours, QPO frequency and X-ray lu- 
minosity are in general correlation with each other, sources that 
span more than two orders of magnitude in luminosity show 
kHz QPOs that cover the same range of frequencies. The frac- 
tional rms amplitudes of the kHz QPOs are ^2 — 20 times 
larger in the atoll sources, at luminosities of ~ 0.005 to ~ 0.2 
of the Eddington luminosity, than in the Z sources, which are 
near-Eddington, and hence 5 to 400 times more luminous. For 
instance, the rms amplitudes of the lower and upper kHz QPO 
in the Z source GX 5-1, which is thought to emit at the Edding- 
ton luminosity, are 1 — 3 % rms and 2 — 4 % rms, respectively, 
whereas for the atoll source 4U 0614+09, with a luminosity be- 
low ~ 0.01 of the Eddington luminosity, the amplitudes of the 
lower and upper kHz QPOs are 6 — 25 %, and 6 — 18 %, re- 
spectively (Wijnands et al. 1998b; van Straaten et al. 2000, 
see Hasinger & van der Klis 1989 for a definition of the Z and 
atoll classes). In most cases, the decrease of the fractional rms 
amplitude of the kHz QPOs in the Z sources compared to the 
atoll sources is less, often by large factors, than the correspond- 
ing increase in luminosity which, as in the case of the paral- 
lel tracks within individual sources, once again argues against 
the idea that the difference in luminosity at constant QPO fre- 
quency is only due to changes in some luminosity component 
that does not affect the QPO. However, while there might be a 
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parameter that changes from one source to another to compen- 
sate for the dramatic increase in luminosity leaving the QPO 
frequency unchanged (e.g., the magnetic field strength; White 
& Zhang 1997), this can not be the explanation for the discrep- 
ancy within individual sources discussed in this paper. 

In principle, our results could be explained if the kHz QPO 
modulates the luminosity from the radial inflow. However, it 
is unlikely that this is the case. Lee & Miller (1998) have re- 
cently suggested that it might be possible for the Comptonizing 
corona to change the rms amplitude of the kHz QPOs indepen- 
dently of the modulation of the radiation emitted by the disk 
(see also Miller et al. 1998; Lehr et al. 2000). If this were the 
case, changes in the properties of the corona in different tracks 
of the v — ix diagram might be the reason for the insufficiently 
rapid decrease of rms amplitude of the QPOs with count rate 
(Fig. 6). (Although it is irrelevant to the present issue, we note 
in passing that some of the conclusions in the paper by Lee & 
Miller are based on the wrong sign of the QPO phase lags; see 
Vaughan et al. 1997, 1998). 

However, it is known that in 4U 1608-52 (and in several 
other sources) at fixed QPO frequency, the hard color does not 
change significantly between different tracks in the v — Ix dia- 
gram (Kaaret et al. 1998; Mendez et al. 1999; Mendez 2000). 
On the other hand, here we have shown that the rms spectrum 
of the lower kHz QPO in 4U 1608-52 (and of the only kHz 
QPO in Aql X-l) does not depend upon the track that the source 
occupies in Figure 2. As long as the high-energy part of the X- 
ray spectrum is produced by inverse Compton scattering in this 
corona, and the corona is also responsible for the rms spectrum 
of the kHz QPO, these two results indicate that the properties 
of the corona do not change significantly from one track to an- 
other. This argues against the idea that changes in the prop- 
erties of the corona could explain the discrepancy between the 
observed and the expected rate at which the rms amplitude of 
the kHz QPOs decrease with count rate (Fig. 6). 

Our simple scenario does not rule out any of the QPO models 
so far proposed, for none of these models explicitly predicted 
the parallel track phenomenon or the behavior of the QPO rms 
amplitude across the parallel tracks. But given that several of 
the QPO models are compatible with the simple scenario de- 
scribed in this paper, our results rejecting this scenario can be 
used to further constrain those QPO models. 

Recently, a new mechanism has been proposed that could ex- 
plain the parallel tracks in the v — Ix diagram of single sources 



as well as across the ensemble of sources. Instead of two physi- 
cal parameters (e.g., M D and Mr) that vary independently, van 
der Klis (2001) proposes that source luminosity and kHz QPO 
frequency are both determined by a single parameter, to which 
the system has both an instantaneous and a long-term average 
response, van der Klis (2001) discusses scenarios where this 
single parameter is the disk accretion rate Md, and QPO fre- 
quency (set by the inner disk radius) is determined by the ratio 
of Md to the total luminosity. He shows that he can repro- 
duce the observed parallel tracks and several of their observed 
characteristics if the long-term average response to Md is, for 
instance, Mr, and luminosity depends upon the total mass ac- 
cretion rate. In models of this type it turns out that QPO fre- 
quency depends upon Md normalized by its own long-term av- 
erage (in our example, Mr). So, points with the same QPO 
frequency on different tracks correspond, in this scenario, to 
source states where both Md and Mr are different, but their 
ratio is the same. The relative strengths of the QPOs in the two 
cases will depend on the strengths with which they are formed 
(at the same disk radius but under conditions of different ac- 
cretion rate), and on the propagation effects of the QPO signal 
through radial flows of different density. This model predicts 
no simple linear dependence on the amount of "extra" X-rays. 
If formation strengths are similar then a higher-density radial 
flow, which (dependent on its speed) would probably occur at 
higher luminosity, might suppress the QPOs more. 
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